Numerical renormalization group study of the correlation functions of the 
antiferromagnetic spin-i Heisenberg chain 
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We use the density-matrix renormalization group technique developed by White jl] to calculate 
the spin correlation functions < S^^iS^ >= {—1)''lo{1,N) for isotropic Heisenberg rings up to 
A*' = 70 sites. The correlation functions for large I and A'^ are found to obey the scaling relation 
Lu{l,N) = Lo{l,oo)fxY{^/^) proposed by Kaplan et al. , which is used to determine uj{l,oo). The 
asymptotic correlation function tj(Z, oo) and the magnetic structure factor S{q = n) show logarithmic 
corrections consistent with uj{l, oo) ~ aV^ncl/l, where c is related to the cut-off dependent coupling 
constant 5e//(^o) ~ l/ln(c/o), as predicted by field theoretical treatments. 



Although the exact ground state of the spin- 1/2 chain 
is explicitly known, the Bethe-Ansatz wave function is 
far too complex to derive directly the spin-correlation 
functions. Other methods like bosonization or conformal 
field theory have to be used to get information about 
the asymptotic behavior of these functions. In general 
quantum spin chains and in particular their continuum 
versions are very active fields of research, because they 
serve as a testing ground of various analytical approaches 

. Recent field-theoretical studies predict the existence 
of logarithmic corrections to the finite-size scaling of the 
energies of these systems and also to the power law be- 
haviour of the spin correlation functions stemming from 
marginally irrelevant operators . The modification of the 
power law originally derived by Luther and Peschel 
has its physical origin in Umklapp scattering processes 
which appear in the fermionic representation of the model 
after Jordan- Wigner transformation and has been antici- 
pated a few years ago . Logarithmic corrections to the 
scaling of the energies and to the correlation functions 
were obtained by Affleck et al. |6) applying conformal 
field theory to Wess-Zumino-Witten non-linear-u models. 
Giamarchi and Schulz |^] and Singh et al Q used a renor- 
malization technique to study the sine-Gordon Hamilto- 
nian and obtained to leading order (— l)'(ln for 
the asymptotic decay of the spin correlation function in 
the case of the isotropic Heisenberg chain. 

In spite of this analytical progress, numerical attempts 
have given contradictory results. Kubo, Kaplan and Bo- 
rysowicz M found a small logarithmic correction of the 



form (Inl)'^ /I with an exponent 0.2 < a < 0.3 instead of 
0.5 as predicted by theory. Later attempts to check the 
theoretical prediction were made by Liang [ p^ and by 
Lin and Campbell , who reported the absence of log- 
arithmic corrections (cr = 0) for spin-1/2 chains, whereas 
Sandvik and Scalapino | [T^ report an exponent a ~ 1/2. 

Progress on the numerical side was hampered mainly 
because highly accurate diagonalization results could be 
obtained only up to iV = 30, while Monte Carlo data for 
larger systems had too large statistical errors. An equally 
important reason was that it had been assumed that the 
data can be analyzed using the universal asymptotic law 
(Iniy /I P-p^. In this paper we clarify and conclude this 
long dispute on the numerical evidence of the logarithmic 
corrections by considering the nonuniversal scaling of the 
coupling constant. 

We use the density matrix algorithm (DMA) |l| to 
study the large-distance decay of the correlation func- 
tions and find that they can be calculated with such high 
precision for sufficiently large systems that the subtle log- 
arithmic corrections to the correlation functions can be 
resolved. This technique leads to highly accurate results 
for much larger systems than those which can be solved 
by straightforward exact diagonalization. The DMA al- 
lows for a systematic truncation of the Hilbert space by 
keeping the most relevant states in describing a state 
{e.g. the ground state) of a larger system, instead of the 
lowest energy states usualy kept in previous real space 
renormalization techniques. A general iteration of the 
method consists of: i)The effective Hamiltonian defined 
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for the superblock l+2+r+2' (where the blocks 1 and 
1' come from previous iterations and blocks 2 and 2' are 
new added ones) is diagonalized to obtain the ground 
state \ip > (other states could be also kept), ii) The 
density matrix pai = ipij^'i'j is constructed, where 
i'i] — < i® >, the states \i > {\j >) belonging to the 
Hilbert space of blocks 1 and 2 (1' and 2'). The eigen- 
states of p with the highest eigenvalues (equivalent to the 
most probable states of blocks 1+2 in the ground state of 
the superblock) are kept up to a certain cutoff, iii) These 
states form a new reduced basis to which all the opera- 
tors have to be changed and the block 1+2 is renamed as 
block 1 . iv) A new block 2 is added (one site in our case) 
and the new superblock (l+2+l'+2') is formed as the 
direct product of the states of all the blocks (the blocks 
1' and 2' are identical to blocks 1 and 2 respectively). 

The method has been applied successfully to several 
problems such as the Haldane gap of spin-1 chains 
the one-dimensional Kondo-insulator ||l^ and the two- 
chain Hubbard model [T5| . 

We used the DMA method keeping up to 200 states per 
block, one target state (the ground state) and periodic 
boundary conditions. The ground state energy has a rel- 
ative error of ^ 10^^ for a system with = 70 sites and 
~ 10~^ for N = 100, as compared to the exact finite-size 
energies calculated using the Bethe-Ansatz. The error 
in the spin correlation function for large distances is al- 
most two orders of magnitude larger than the error in 
the energy, as we estimated by comparing our calcula- 
tions with exact results for 30 sites |l^ . Because of this 
we present results up to A^ = 70 so as to have a large 
enough accuracy. 

In the following we consider the average value 
uj{l, N) = -1,N) + 2uj{l, N) + uj{l + 1, A^)] (where 

< ^n+i^n to remove the even-odd-Z 

oscillations of the spin correlation function, which are of 
order in ui and ^ in w To extract the val- 

ues of the correlation function for the infinite system, we 
adopt the scaling relation of Kaplan et al. |0] 



id{l,N)=uj(l,^)f{l/N), 



(1) 



which is expected to hold for sufficiently large I and A^. 
The scaling function is given hy f(l/N) = fxY^il/N) in 
terms of the scaling function of the XY-model 

fxvix) = 1 + .28822 sinh2(1.673a;). (2) 

To check the scaling we plot in Fig. 1 

Z{l,N)::^uj{l,N)/u;{l,2l) (3) 

as a function of l/N. We find that all the curves coincide 
within 0.3% in a unique curve for I > 7, which confirms 
the scaling hypothesis. Following the proposal of Ref. |^] 
we performed a least square fit for these curves using the 
function 



Z{l,N) = 



fxY{l/N) 



fxYil/2) 



(4) 



and obtained for the exponent a — 1.8052 (the small 
error is determined from the dispersion of the fitted ex- 
ponents considering all the curves). The fit is also shown 
in Fig. 1. This implies that the scaling function defined 
in Eq. 1 is: 



fix) 



1 



.28822 sinh^(1.673a;)] 



1.805 



(5) 



The surprising quality of the scaling is shown in Fig. 2 
where we plot the averaged N) as a function of l/N 
for several values of I. From Eq.l we see that the size 
dependence of the correlation function is given by the 
function f{l/N) rather than the normally used 1/A^^ 
behaviour. We have fitted these curves with a one- 
parameter fit using Eq. |^ with (D(?,oo) as the free pa- 
rameter. We use these extrapolated values a)e(Z,oo) to 
study the logarithmic corrections below. 

There are very small deviations from scaling which can- 
not be seen in Fig. 2 , which look systematic and are prob- 
ably not due to the Hilbert-space truncation. They ap- 
pear as a small A^-dependence of ui{l, N)/ f{l/N). These 
deviations are less than 0.4% of the value at l/N —1/2. 
We note that these deviations could induce an error on 
uj{l, oo) of at most 1% for the largest /-values. 

To display the logarithmic behaviour of the correlation 
function we have plotted lu:{l,N)/f{l/N) as a function 
of Z (/ > 7) for different A" values in Fig. 3. Here we 
can see the small A^-dependence for a given I mentioned 
above, but a clear logarithmic behaviour is found. We 
also show the extrapolated values luJe{l, oo). To visualize 
the magnitude of the finite-size correction to the corre- 
lation functions, we show in the inset the bare data for 
a finite system, loj{l,N) as compared to the finite-size 
corrected values loj{l, N) / f{l/N). We see that the cor- 
rections are largest for I ~ A^/2 (as can also be seen from 
Eq.|) 

The logarithmic corrections to the correlation function 
follow from the scaling relation 



G,(r) ro 
__ = _ exp 

G^{ro) r 



-Anb, / d{lnr')g{r') 

Jro 



(6) 



where we use the notation of Ref. 6 with Gi (r) = Lu{r, oo) 
and / replaced by the continuous variable r. The r- 
dependence of the coupling constant to one-loop order 
is given as 



giro) 



1 + nbg{rn) ln(r/ro) ' 



(7) 



For spin- 1/2 the parameters are determined by 6 = 4/\/3 
and Abi/b = —1/2. The r-dependence of the coupling 
constant leads after integration to the multiplicative log- 
arithmic correction of the form 
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Gi{r) 
G.(ro) 



9{r) 



(8) 



For the following analysis we insert Eq.(7) and obtain 



ln(- 



ro 



(9) 



with 5e//(^o) = 47r(7(ro)/\/3- This equation shows that 
the universal, i.e. coupling constant independent, asymp- 
totic relation ^(Inr)^/^ is only reached for sufhcicntly 
large distances r. 

We fitted the curves in Fig. 3 with the dependence 
predicted by renormalization group 



Ilu{1, oo) = a-\/ln(cZ), 



(10) 



exp(l/(7e//(^o))/'o defines the scale on which 

VW)/^ is 



where c 

the asymptotic behavior uj{l,oo) ~ y'ln[t)/l is ap- 
proached. Figure 4 gives a comparison of the data and 
our best fit with a ^ 6.789 x 10"^ ^nd c = 23.21 . This 
implies that the asymptotic regime, i.e. ln(c)/ln(/) << 
1 , is reached only for chains with more than several thou- 
sand sites. The small deviations may either stem from re- 
maining uncertainities in the determination of <De(^, oo) or 
from the one-loop calculation, which is exact only to or- 
der 0{g^). From the value for c we obtain for the coupHng 
constant ge//(20) = 0.163 at Iq = 20. This may be com- 
pared with an effective coupling constant g^ff {20) ~ 0.26 
deduced numerically from the scaling of the ground state 
energy and triplet and singlet excitations 

We stress that our analysis of the data differs from 
earlier numerical studies where the asymptotic expres- 
sion <D(Z,oo) = a(ln^)°'/Z was considered assuming c = 1 
and with cr as a free parameter [p|-|ll|. Sandvik and 
Scalapino on the other hand suggested that previ- 
ous numerical studies did not succeed in finding 
the proper exponent since the scaling relation (1) may 
not hold in presence of logarithmic corrections. They 
proposed instead an alternative relation which connects 
uj{l, N) with the asymptotic correlation function uj (I, oo) 
and which does not obey (1) [|l^. After subtracting the 
oscillatory 1/^^-contribution from the correlation func- 
tion they analyse the ratio D{1, N) = uj{l, N)l/\/\n I, as- 
suming the log-corrections are given by VhiT. Given our 
results D{l,N) cx f{l/N)^/h\ cZ/\/ln I, i.e. the functional 
used in Ref.l2 accounts for finite size effects but also for 
parts of the log-corrections, consequently their analysis 
concerning the form of these corrections is not conclusive. 

Finally it should also be possible to determine the 
logarithmic corrections from the TV-dependence of the 
structure factor S'Ar(q = tt) = ■ Earlier at- 



tempts [P Hliyi8| were not successful because this expres- 
sion also involves large finite size effects as uj{l,N) ^ 
u){l^oo)f{l/N). It is therefore better to consider the 
quantity 



S'7v(7r) 



N/2 

E 

-N/2+1 



^{l,N) 
fil/N)- 



(11) 



Given Eq. [l^ for oo) at large / one expects 57v(7r) = 
const + ^alTi^^^{cN/2). Fr om our fit of Sat wc find val- 
ues for a and c which are very close to the parameters 
deduced from oo) (see Fig. 5). 

In conclusion we have shown that the scaling relation 
proposed in |^ provides a very accurate description of 
finite-size effects for large enough systems and distances 
(l > 7), and we have obtained an improved value for the 
exponent entering the expression for the scaling function. 
Furthermore we have shown that the correlation function 
uj{l, oo) of the infinite system, which can be determined 
from A^-site rings for I < N/2, does not obey the uni- 
versal asymptotic law (ln^)°'/^ as was assumed in previ- 
ous numerical work, but is governed by the nonuniversal 
scaling of the coupling constant. Our data confirms the 
multiplicative logarithmic corrections to the spin correla- 
tion function of the form uj{l, oo) ~ a\/ln cl/l as derived 
from quantum field theory [^-||, and we determine in 
particular the scale parameter c of the logarithmic term 
and the related effective coupling constant which has not 
been obtained by these approaches. 
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FIG. 1. The ratio Z{l,N) = uj{l,N)/uj{l,2l) versus l/N 
for Z > 7 and 14 < iV < 70. The scaUng function (4) with 
a = 1.805 is shown as solid line. 

FIG. 2. Correlation function lij{1,N) versus l/N for 
Z = 7 — 29. The solid curves are single parameter fits us- 
ing the scaling relation (1) with Q{1, oo) as free parameter. 



FIG. 3. Comparison of lu!{l,N)/f{l/N) versus / for 

A'" — 14 — 70 with toe(',oo) (solid curve). In the inset 
lid{l,N)/f{l/N) (filled circles) and lu!{l,N) (open circles) are 
shown for N = 70 showing how f{l/N) corrects the finite-size 
effects. 

FIG. 4. Logarithmic corrections are clearly seen in the 
quantity lLUe{l,oo) versus /, which is well fitted by the an- 
alytical expression (10) (solid curve). 

FIG. 5. Sjv(7r) (Eq. 11) vs. ln{N). Comparison of data 
(+) and fit (solid line) using Sjv(7r) = const + |aln^/^(cAr/2) 
with const = -4.06 10-^ a = 6.67 10"^ and c = 25.5. 
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